ASYMPTOTICS FOR SOLUTIONS OF ELLIPTIC EQUATIONS IN DOUBLE 

DIVERGENCE FORM 

VLADIMIR MAZ'YA AND ROBERT MCOWEN 

Abstract. We consider weak solutions of an elliptic equation of the form didi(aiju) = and 
their asymptotic properties at an interior point. We assume that the coefficients are bounded, 
measurable, complex-valued functions that stabilize as x — » in that the norm of the matrix 
(aij(x) — on the annulus B2 r \B r is bounded by a function Cl(r), where £2 2 (r) satisfies 
the Dini condition at r = 0, as well as some technical monotonicity conditions; under these 
assumptions, solutions need not be continuous. Our main result is an explicit formula for 
the leading asymptotic term for solutions with at most a mild singularity at x = 0. As a 
consequence, we obtain upper and lower estimates for the L p -norm of solutions, as well as 
necessary and sufficient conditions for solutions to be bounded or tend to zero in L p -mean as 
r -» 0. 



0. Introduction 

We are interested in the local behavior of weak solutions to the elliptic equation in "double 
divergence form" 

(1) Au := didj(aij(x)u(x)) = 0, 

where we have used di = d/dxi and the summation convention; the coefficients Qtij — CLj % clX6 
bounded, measurable, complex-valued functions in a domain to be specified. The operator A 
arises naturally as the formal adjoint C* of the operator in "non-divergent form," 

(2) C = a,ij(x)didj. 

Solutions of are not only important for the solvability of Lu = /, but for properties of the 
Green's function for C. When the coefficients dy are real-valued functions, the operators C and 
A have been studied by Sjogren Bauman (|3], g], Fabes and Stroock ([T2]), Fabes, 

Garofalo, Marin-Malave, and Salsa Escauriaza and Kenig (^HDi ancl Escauriaza ([5], [5]); 

these papers use techniques, such as the maximum principle, that rely on the coefficients being 
real-valued. Not only do our techniques apply to complex-valued coefficients, but they give 
additional information for real- valued coefficients; we shall explain this in some detail at the end 
of this Introduction. 

We want to study weak solutions of in a neighborhood of an interior point of the domain, 
say x = 0, where the coefficients ay approach <5y in the sense that 

(3) sup |a(x)-I| <0(r); 

r<|x|<2r 

here sl(x) is the matrix (ay), I is the identity matrix, | • | denotes the matrix norm, and Q(r) — * 
as r — * in a manner that we shall describe. We remark that, when the coefficients are real- 
valued, the more general case obtained by replacing 5ij by constants ay that satisfy the ellipticity 
condition can be reduced to © by means of an affine change of the x variables. Of course, this 
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reduction of the more general case to J2J is not available when the constants oiij are complex- 
valued, but we have chosen to treat the special case otij = 8%j in order to take advantage of 
technical simplifications in the formulations and proofs of our results. 

The specific hypotheses that we impose on the function f2(r) in are as follows: 

(4) jf 1 ^l d t < co, 

(5) Q( r ) r~ 1+£ is nonincreasing for < r < 1 and, 



(6) r ™ 6 is nondecreasing for < r < 1; 

here e > 0. Clearly, 10} together with JSJ) or Jfj} implies that O(r) — > as r — > 0, so the coefficients 
Ojj are approaching 5y as a; — > 0, although perhaps at a slow rate. 

A weak solution of QJ in a domain U C R™ is a function u G L\ oc (U) that satisfies 

(7) / dyO) u{x) djdifjix) dx = for all r? G C£°(E/). 

Weak solutions of Q need not be continuous under our assumptions on the coefficients, so to 
measure growth or decay as x — > 0, we will use the mean in L p for some p£ (1, oo): 

(8) M p (w,r) := (P ^ \w\ p dx^j '\ 

where A r is the annulus B2 r \B r with B r = {x G K : \x\ < r}; here (and elsewhere in this paper) 
the slashed integral denotes the mean value. (We will also use the notation M p (w, r) when w is 
vector or matrix valued; in this case, \w\ denotes the norm of w.) 

For our local results, we will consider in the unit ball B\ and we will assume 

Jo 1 

where 5 is sufficiently small. In fact, this represents no additional assumption on f2(r) since we 
could replace B\ in what follows by a very small ball B 1 in order to make the integral tt 2 (t)t~ 1 dt 
as small as necessary. 

At times it will be useful to consider solutions of in all of W 1 ; in that case, we assume that 
a ij = ^ij outside of B\ . Our first result concerns such a solution. 

Theorem 1. Let n > 2, p £ (l,oo), and tt(r) satisfy |3|) ; and J^. There exists a weak 
solution Z G Lf oc (M. n ) of equation in R™ satisfying 



(10) Z(x) = exp 

where £ satisfies 



2x d V 



/ {au{y)-na l] {y)y l y j \y\ ) 

\oBi\ J Bl \B M \y 



M p (C,r) < c max (fl(r), J dtj for0<r<l. 

Remark 1. T/ie proof of Theorem 1 shows that Z has a limit at infinity, Z{po), satisfying 

\Z{x)-Z(ao)\ <cVS\x\- n . 

But we are more interested in the behavior of Z at the origin; Corollary 3 below can be used to 
show that Z has at most a mild singularity at the origin (cf. \T^j ). 
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Our second theorem uses the solution Z from Theorem 1 to characterize the asymptotics (as 
x — > 0) of weak solutions of 0); because this is a local result, we consider a solution in B±. 

Theorem 2. Let n > 2, p S (l,oo), and £l(r) satisfy J3|), and 0). Suppose that u S 

^ 7 ioc(-^ 1 \{^}) * s a wea ^ solution of in £?i subject to the growth condition 

(11) M p (u,r) < cr 2 - n+eo , 

where Eq > 0. TTien £/iere exists a constant C (depending on u) such that 

(12) u(a;) = CZ(x) + «;(ar), 
where the remainder term w satisfies 

(13) M p (w,r) < cr 1 ' 61 
for < r < 1 and any Ex > 0. 

Remark 2. TTie restriction n > 2 in Theorem 2 is caused by the existence of solutions for the 
Laplacian with logarithmic growth at x = when n = 2. A refinement of the techniques used in 
proving Theorem 2 would be required to cover the case n = 2. 

The following two results are immediate consequences of Theorems 1 and 2. 

Corollary 1. Under the hypotheses of Theorem 2, the condition 

(14) lirninf / Re(a u (y) - na ij (y)y i y j \y\~ 2 )-- - > -co 

r > j B\\B r \y\ 

is necessary and sufficient for the boundedness of M p (u,r) as r — ► for all u € Lf oc (Bi\{0}) 
which are weak solutions of and satisfy 

Corollary 2. Under the hypotheses of Theorem 2, the condition 

f ^y 

(15) lim / Re(au(y) - na tJ (y)y i y j \y\- 2 )—- = +oo 
r ^°JB 1 \B r \y\ 

is necessary and sufficient for M p (u,r) — ► as r — > /or a// u 6 Lf oc (Bi\{0}) which are weak 
solutions of and satisfy 

Moreover, we can use Theorems 1 and 2 to derive upper and lower estimates for the L p -mean of 
solutions. Since (Re(ajj) — <%) is a symmetric real matrix, it is not difficult to verify that 

(16) _ 2 (n _ 1)| Rea(n) - I| < Re(a u (y) - na lJ (y)y l y J \y\- 2 ) < 2(n - 1)| Rea(y) - 1 1. 

Using | Re a(y) — 1 1 < | a(y) — 1 1 and (JSJ) , we find that IjlOl) also yields the following upper and 
lower estimates. 

Corollary 3. The solution Z in Theorem 1 satisfies for r £ (0,1) 

(17) ciexp (-2(n- 1) J fl(s)yj < M P (Z, r) < c 2 cxp ^2(n - 1) J 0( s )y 
where C\, c 2 are positive constants. 

Even when the coefficients a^- are real, the upper and lower bounds l|17|l appear to be new. 

The principal analytic content of our results is contained in Theorem 1. The method of its 
proof is independent of, but related to, the asymptotic theory developed in ^7]. In particular, 
Lp-means of type © were extensively used in ^5] and ^B]. The asymptotic formula that we 
obtain is analogous to that of |18j . where an asymptotic representation near the boundary was 
obtained for solutions to the Dirichlet problem for elliptic equations in divergence form with 
discontinuous coefficients. 
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Now let us turn to the comparison of our results with the extensive work of the authors cited 
in the first paragraph; we refer to the excellent exposition in Escauriaza [H] for a more detailed 
description of these results and references to the literature. When C is uniformly elliptic with 
real- valued, measurable, and bounded (although not necessarily continuous) coefficients on R™, 
the previous work shows: i) any nonnegative weak solution Z S L\ oc {B{) of QJ in B\ satisfies 
Z e L n '^ n ~ 1 \Bi) (cf. 12 ); ii) for a fixed nontrivial nonnegative weak solution Z 6 L\ oc (B\) 
of in Si, and every weak solution u G L\ oc (B{) of in Si, the function u/Z is Holder 
continuous in B\ (cf. [3]); and iii) the existence of a unique nonnegative weak solution Z of Q 
in R™ that satisfies J B Z dx = \Bi\ and may be used to estimate the Green's function for the 
operator C on R n (cf. |§]). These previous results are quite general, but do not apply to the case of 
complex coefficients that we consider because they rely upon the maximum principle. Moreover, 
even when the coefficients are real-valued, our results are somewhat different in nature than the 
previous ones: instead of estimates, we have obtained an asymptotic description of Z near a point 
where the <Zjj are continuous in the sense of (|3J). For a more direct comparison, when the are 
real-valued and continuous at 0, Escauriaza |5] has obtained upper and lower estimates for the 
L 1 -norm of Z: for any 7 > there exists a constant iV 7 such that 

(18) N~ 1 r'<j' Zdx<N~ l r-~ l for < r < 1. 

When the ay satisfy our stronger sense of continuity the explicit upper and lower bounds 
<|17fl imply L p -bounds analogous to l|18l) . In fact, as we shall see in Section 2, the assumptions on 
f2(r) imply that f2(r) < cs/S for < r < 1 (cf. I@B}), so we easily obtain from (|17|) 

(19) < M p (Z,r) < iV 7 r~ 7 for < r < 1. 

More importantly, however, our necessary and sufficient conditions for solutions to be bounded 
in L p -mean as r — > (Corollary 1) or tend to zero in L p -mean as r — > (Corollary 2) have not 
been obtained previously. 

1. Preliminary Estimates 

In this and the next section, we will use the spherical mean of a function w. For notational 
convenience, we denote the spherical mean using an "overbar" : 

(20) w(r) = f w(r6)ds. 

J dBi 

This should cause no confusion with complex conjugation since we will not have occasion to use 
the latter in these two sections. In particular, in this section we are concerned with solving an 
equation of the form 

(21) -Av = did^) - didjiFij) mW 1 . 

Here Fij £ Ll oc (M. n ) and derivatives are interpreted in the sense of distributions. The norm of 
the matrix T = {Fij) will be denoted by 

Proposition 1. Suppose that F^ £ L p oc (R™\{0}) satisfies 

(22) f \F{x)\dx+ I |^(a;)||a;r n - :L (ix<oo. 

J\x\<l J\x\>l 

Then there exists a weak solution v € L p oc (R™\{0}) of }21\) that satisfies 

(23) M p (v,r) <c\M p (F,r)+r f \F(x)\\x\- n - x dx + r~ n [ \F{x)\dx ] , 

\ J\x\>r J\x\<r J 
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where c is independent of r and we have introduced 

i/p 



M p {w,r) := I f \w(x)\ p dx 

\J r/2<\x\<4r 



r/2<\x\<ir 

Proof: It suffices to prove the result for Fij g C^°(R"\{0}) since the general case can be handled 
by a standard approximation argument. The function v is dchned by convolution with V, the 
fundamental solution for the Laplacian: 



v = T * (didjFij - didjFij). 
Using f Rn f{y)g(\y\)dy = J* R „ J(\y\)g(y) dy, we can write this as 

v(x)= f (T(\x - y\) - T{\x - -|)(|2/!)) didjF^y) dy. 

Now to compute the spherical mean of the fundamental solution, we can use the mean value 
theorem for harmonic functions to conclude that 



r(|^.|)(M) = r(max{M,M}) 

This enables us to express v as 

d 2 f .. , , . „„ „ f i)-'F, , r „„ „ -•'-/•;, 



Now, integration by parts yields 



and 

}2 



v{x) = oxh L r(| * - ^ (y) dy - r(N) L., IS dy - L >,., r(M) S * 

/ ^dy=( ^l^dS 
J\v\<\x\ dyid Vj J\ y \=\ x \ d Vj \v\ V ' 

I T (\y\)^ L dy = - [ T'{\y\)r-^dy- I T{\y\)^l^-dS y 
J\v\>\x\ hv\>\A \V\ Q Vj J\v\=\*\ ^ W 

= / 7T (r'(\y\)Tl) F iMdy+ I (v\\y\)^F ii {y)-T{\y\)^)^-dS y . 
J\v\>\x\dyj\ \y\J J\v\=\x\\ \y\ dy 3 J\y\ 

Consequently, 

-T'{\x\) I y^i Fij (y)dS y . 

J\v\=\x\ \y\ 

Now introduce xo and Xoo as the characteristic functions of B r / 2 and S| r , and let xi 
1 — xo — Xoo be the characteristic function of the annulus A r := B4 r \B r / 2 - Then 

v{x) d^~dx~ 3 L T{lx ~ yi){xiF * M dy = 

I il^r^y^^Mdy-T'^) f y -^±F 3 {y)dS y . 
Jb 4 ,.\b 1x1 \oyioyj J J\y\=\x\ \y\ 
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We can estimate the four integral kernels and obtain that the right hand side is bounded by 

WiM\dy 



This provides us with the following pointwise bound 

d 2 



n+1 



v(x) 



dxidxj 



T{\x-y\){ X iFi 3 ){y)dy 



(24) 



dy 



<cl|^|(N) + N J^\F ij {y)\-^ T + \x\- n J B \F l3 (y)\dy 



Using the ZAboundedness of singular integral operators on 1™ (see |27|h we have 
d 2 



dxidx 



T(\x - yDixiF^dy 



LP(A r ) 



< ||X1^||zp(R») = \\HLv{A r y 



Elementary estimates may be applied to the remaining terms in (|24[) to obtain J22J), completing 
the proof. □ 

The integrals in (|23Jl can be estimated in terms of M p . For example, we substitute 



bl-"" 1 =c„ 



'-<\y\<\x\ 

into the first integral and change order of integration to obtain the estimate 

rf ^{xm-n-Hx < cr [ [ \f {x )\ dx ^L^< cr f°° Mp {F,p)%. 

J\x\>r J\y\>r/2J\y\<\x\<2\y\ \V\ ^ Jr/2 P 

Similarly, we can show 

r~ n [ \T{x)\dx<cr~ n ( M p (T, p) p n ~ l dp 

J\x\<r JO 



and 



4r 



M P {T, r)P < cr- n / M P [T, p) p p^ 1 dp. 

Jr/2 

Elementary estimates show that terms involving integration over r/2 < p < r and 2r < p < 4r 
can be respectively dominated by the terms involving integration over < p < r and p > r, so 
we obtain the following. 

Corollary 4. Under the hypotheses of Proposition^ the weak solution v obtained there satisfies 



(25) M p (v, r)<c[r j M p {F, p)p- z dp + r~ n J M p {F, p)p n - i dp 
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2. Proof of Theorem 1 



We shall prove Theorem 1 by reducing the problem of finding Z to solving an operator equation 
of the form (I + T)V = f, where V and / are elements of a Banach space X of functions on 
M n \{0}), and T is an integral operator of small norm on X. However, this reduction will take a 
few steps. To begin, let r = \x\, 9 = x/\x\, and r\ € Qj°((0,oo)) be arbitrary. For Z € Lf oc (M. n ) 
to be a weak solution of QJ, we must have 







V "(r) / Z(r0) aij (r6 



%9 jds 



didj7](\x\)a,ij (x)Z(x) dx 
r)'(r) 



Z(r6)(au(r6) - ay(r6 



iob 1 r 
where ds denotes surface measure on the unit sphere, 8B%. Hence, 



ds ) r n ^dr. 



= 



Z{rff)ais{rff)Qi6jd8 



Z{r6){a H {r6) - a lJ (r9)e i 9 J )ds dr, 
asi / 

where the derivative is understood in the distributional sense. This implies 



(26) 



r n-l± 
dr 



Z(r9)a ij (r9)9 i 9 J ds + r ri 



OB! 



Z(r9)(au{r9) - na^ ^9)9,9^8 = C, 



where C is an arbitrary constant. In what follows, we will take C — 0; as we shall see, the solution 
that we construct will in fact be a weak solution of Q on all of K™, not just M™\{0}. (See also 
the Remark at the end of this section.) 
Let us introduce 

(27) v{r9) := Z(r9) -Z(r), 

where Z is the spherical mean as in H20|) . We may now express (12611 as 

(28) 



y>M + ^-y{r) = -Kv(r), 
r r 



where 
(29) 

and 
(30) 
with 
(31) 



y{r):=y Z{r6)a ij {re)e i e j ds, 



dBi 



Q(r) :=n — 



in (|28() we also have used 
(32) Kv{r) := 



a(r) 

n (l ( ) := J (i ri {r9) ds am i a{r) : 



dBx 



%6j ds; 



v(r9)au(r9) ds 



a(r) 



v(r9)a l] {r9)9 t 9 j ds. 



It follows from © that |ctoM — n \ < c ^( r )i \ a ( r ) — 1| < c ^( r )j an d 
(33) |Q(r)| < cfl(r), 
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so Q(r) — > as r — > 0. Since v(r) — 0, we can also write (|3(J|) as 

Kv{r) = !f v(r6)(a u (r0) ~n)ds- v(r9)(a l3 (r9) - 6, i )D,0 j ds. 

J OB! a V) J 9Bi 

In this last form it is evident that K satisfies 

(34) M p (Kv,r) < C n(r)M p (v,r). 

To obtain another equation involving y and v, we start from the identity 



(35) Av = didj((aij - 5ij)v) - didj{{a io - - 5 i3 )v) + didj((a,ij - 5 l3 )Z) - didj((a,ij - 8 lJ )Z). 
Noting that 



y{r) = a(r)Z(r)+f v{r6)(a l3 (r0)Mj - I) da, 

dBi 



we can rewrite (|3*51) as 



(36) Av = didj{Bij{v)) - didj(B i3 (v)) + didj^y) - didj^y), 
where 

(37) B l3 {v){x) := ( aij (x) ~ <%) (v{x) - -^-rj' u(r0)(ay (rO)^ - 1) ,/* 
and 

(38) fajix) 



a(r) 

Using 0, it is clear that 

(39) M p (£? 4 »,r) < cQ(r)M p (v,r) 
and 

(40) |<^0)| < cO(r). 

These estimates indeed hold for < r < oo, where we extend ft to be zero for r > 1. Inverting 
the Laplacian by means of the fundamental solution, (|36[) becomes 

(41) v + Sv + Ty = 0, 
where we may use Corollary 0] with (|39|l to obtain 

(42) M p (Sv, r)<c(r£ n{p)M p (v, p)p- 2 dp + r~ n J Sl(p)M p (v, p)p n ~ 1 dp^j 
and with (|40|l to obtain 

(43) M p (Ty,r) < c (W n(p)M p (y,p)p- 2 dp + r- n J^ Sl(p)M p (y, p)p n - 1 dp^j . 
To simplify the equations, let us introduce the function 

(44) £(r) = expQT^dt), 

where there is no problem with convergence of the integral since Q(t) = for t > 1. Notice that 
E{r) is continuous for r S (0, oo), E{r) = 1 for r > 1, and for any r, p S (0, oo) we have 

(45) J E- 1 (r) J B(p) = exp ( f 9ttl & 
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To derive an estimate for this expression, let us use © with e £ (0, 1/2) and © to conclude 
(46) S> f ^ldt>n 2 (r)r- 2+2s f t x ~ 2e dt > c0 2 (r). 



/2 t J r / 2 



^dt 



< ( - ) for < p < r < 1. 



As a consequence of (|33|) . we have 

(47) (£)^<|exp( ± ^ ( 

Now let us express equations (12811 and (|41(l in terms of the new dependent variables 

(48) Y(r) = E- X {r)y{r) and V{x) = E- X {\x\)v{x). 
Since the operator K only involves integration in 9, it is clear that 

(49) Kv(r) = E(r)KV(r), 
and so the equation f2"%| can be expressed as 

(50) Y'(r) = -KV{r). 

We will be assuming below that M p (V,r) < cQ(r) for < r < 1, so M p (KV,r) < cfi 2 (r), 
enabling us to integrate (|50[l to obtain 

(51) Y{r)=Y(Q)+ f ^Mdp. 



The equation l|41|l . on the other hand, is replaced by 

(52) V + SxV + T{Y = 0, 
where 

(53) Si := E~ X SE and T x := E~ X TE, 

with E representing the multiplication operator defined by the function 

Now let us substitute l|51|l into (|52|l to finally obtain the operator equation that we want to 
solve: 

(54) V + SiV + T 2 V = -F(0)Ti(l), 
where 



E{p) / KV{t)dt 



(55) T 2 V{x) =E- x (\x\)T p _ 

(In notice that T\ operates on functions of a single variable, say p, and Ti(l) represents the 
action of T\ on the function that is identically 1.) For a given choice of Y(0), we can solve 11541 
uniquely for V provided we can show that the two integral operators involved have small norm 
on an appropriate function space. Consider the functions w in L^ oc (M\{0}) for which the norm 

I, I, M p (w,r) M p (w,r) 

(56) w n : = sup f- +sup^= 

0<r<l S2(r) r>l v^r-™ 

is finite, and take the closure to form a Banach space X. We want to show that the right hand 
side of l|541 is in X and that the integral operators Si and T2 map X to itself with small norm. 
It will be useful to observe that the continuity and positivity of |£7(r) | implies that, for any 
w € J Lf oc (R\{0}), we have 

(57) M p (Ew,r) = \E(f)\M p (w,r) for some r = r w € (r, 2r), 
with an analogous statement for E~ x . 
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To show Ti(l) E X, we must estimate M p (Ti(l),r) separately for < r < 1 and r > 1. For 
< r < 1, we can use (15711 to find r £ (r, 2r) so that 

M^T^r) = \E-\r)\M p (T p ^ r [E{p)],r), 

and then use H43[) to estimate 

MpinO), r) < c \E- l {r)\ (r J™ Sl(p)M p (E, p)p- 2 dp + r~ n jf Q(p)M p (E, p)p n - 1 dp) 



= c\E- 1 (r)\\rJ n(p)\E(p)\p-'dp + r- n J n(p)\E(p)\p n -' dp ^ 
where p S (p, 2p) by (|57Jl . But now we can use 1)45(1 and (|47|1 to conclude 

A^(Ti(l),r) < c (r jf il(p)(p/r)°S* p~* dp + r~ n jT Sl{p)(r/p)°^p n - 1 dp) 

Using the monotonicty properties © and 10, we obtain 

M p (Ti(l),r) < c (n(r)r E - cV * £ p- e ~ 1+cVI dp + n{r)r c ^- £ J p^^^dp) . 

Provided S is sufficiently small that e — c\fb > 0, we conclude that 
(58) M p {Tx{l),r) < cfi(r) for < r < f . 

On the other hand, for r > 1 we use i? _1 (r) = f and f2(r) = to estimate 

M p (Ti(l),r) = M p (T^ r [E(p)],r) < cr" n / fi(p)M p (.E, p)^ 1 dp 



= cr- n Q,{p)\E{p)\p n - 1 dp < cr~ n O^p™"^" 1 dp. 
Jo Jo 

Provided n — c we obtain 

(59) Mp(Ti(l),r) < c\/Sr- n for r > i. 

The two estimates (|58|l and <|59[) together confirm that Ti(l) G X. 

Now let us show that Si maps X to itself with small operator norm. We suppose that ||U|| p .q < 
1 and estimate M p (S\V,r) separately for < r < f and r > 1. For < r < 1 we have 
M p (V,r) < f2(r) and we can argue as in the previous paragraph to obtain 

MpiSjV, r) <c(rj tf{p){p/r) c ^ p~ 2 dp + r~ n jT tf(p)(r / pf^ p"- 1 dp) . 

Using H46|l and the monotonicity analysis as above, we conclude that 

(60) M p {S,Vr) £cV - s for0<r<1 

S2(r) 

For r > 1 we use 11470 with r = 1 and l|46[) to obtain 

MptSiV.r) < cr" n / n 2 (p)p"- c ^- 1 dp < cr~ n <5 / p"-^' 1 dp < cr- n 5 
Jo Jo 

provided S is sufficiently small, and we conclude that 

(61) B <cv^ forr>l. 

The estimates (|6U|I and (I61|l together show that Si maps X to itself with small operator norm. 
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Finally, we estimate T%. For < r < 1 and M P (V, r) < fi(r), we argue as before to write 



M p (T 2 V,r) <c(r I n(p) {p/rf^ M p 



KV(t)-,p 



P 2 dp 



P KV(t)^,p 



dp 



But we can use (EU to obtain M p (KV,p) < cO,(p)M p (V, p) < cQ 2 (p) so that 



M p (T 3 V,r)<c(rJ fl(p) [p/r) cV ~ s 



dt 



p 2 dp- 



n(p) (r/p) 



dt 



P n ~ l dp) < c5n(r). 



Using H46|) and the monotonicity argument, we have 

M p (T 2 V ir ) 



(62) 



< cS for < r < 1. 



For r > 1, we use 1)4 7fl with r = 1 and (|33[) to obtain 

M p (T 2 F, r) < cr- n / n 2 (p)p- cVJ +™- 2 dp < cr- n 5 / p -<^+™- 2 dp < cr -n^ 



provided <5 is sufficiently small. These estimates show that T 2 maps X to itself with small operator 
norm. 

Since Si and T 2 have small operator norms on X, we conclude that (|54|) admits a unique 
solution V £ X. Let us now investigate the implications for the weak solution Z{x) — Z(\x\)+v(x) 
that we are trying to construct. Tracing back through the definitions, we see that our solution of 
GJ is given by 

E(r) 



(63) 
where 



Z(x) 



a(r) 



Y(0) + Yi(r) 



Yi(r) 



V(re)a i: j(re)9iej ds + a(r)V{x) , 



KV(p) 



dp. 



Recall that \a(r) — 1| < cfi(r) and V satisfies M p (V, r) < cSl(r) as r — > 0. Moreover, 



M„ 



< c 



o P 



dp, 



so all the terms after the Y(0) inside the parentheses of l|63[) are bounded in M p either by O(r) 
or by J Q r il 2 (p)p~ 1 dp for < r < 1. Let us explore -E(r). Notice that 



Jr P Jr 



-a (p)+na(p)}^+ [' Q(p)[l - a(p)}^- 



P 



\dBi\ j BAB . 
But \Q(p)[l-a(p)]\ < cfl 2 (p), so 



(na ij (y)y i y j \y\ 2 - au(y)) 



dy 



Q{p)[\-a{p)\ 



dp 



f Q{p)[l-a{p)]^ = C- f 

Jr P JO 



Q(p)[l-<*(P)]\-- 
P 



Exponentiating, we obtain 

nl Q(p) 



E{r) = exp 



dp 



C exp 



(l + Co(r)), 
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where |Co(f)| < c /J" Sl 2 (p)p 1 dp. Since we can rescale Z to set CY(0) = 1, we have the formula 
l |10p . This completes the proof of Theorem 1. 

Remark 3. Choosing C — — |c?^Bi | 1 in the proof of Theorem 1 and modifying the rest of the 
argument, we could obtain the asymptotic representation of the fundamental solution of equation 

m 

T A {x) = T(x) + A- 1 d l d J {(a l3 - %)r) + v(x), 

where T is the fundamental solution of A and 

M p (v,r) <c0 2 (r)r 2 "" for re (0,1). 

Here the condition ^ is not necessary; we only need the smallness of Q(r). However, this 
formula, as well as the additional terms in the asymptotic expansion of T_4 can be obtained more 
easily by iteration from the equation Au(x) = 5{x). 

3. Proof of Theorem 2 

Let q = p/{p — 1) and (3,j e R. Let us introduce the weighted L p norm for functions on R" 
with separate weights at the origin and infinity 

(64) \\u\\l q = |K, (Bl) + \\u\\ q LUBC) = [ \u(x)\«\xf«dx + I \u(x)\"\x\^dx, 

' e { lJ j( lJ J\x\<i J\x\>i 

and the weighted Sobolev space W / |'^(R™\{0}) with norm 

(65) J2 \\r lal d a u\\ L ,^ y 

\a\<2 

Notice that L p _ _ 7 (R") is the dual space for L q 0/y (R n ) and the notation WZp P _ 7 (M. n \{0}) will be 
used for the dual of Wp ,q (M. n \{0}) . Many authors have used similar weighted Sobolev spaces to 
study operators like the Laplacian on K™, R™\{0}, and other noncompact manifolds with conical 
or cylindrical ends. 

Using the analysis in [22], [21] or [^Oj, for example, it is easily verified that the bounded 
operator 

(66) A : W£«(R"\{0}) - ^ +2 , 7+2 (M") 

is Fredholm (finite nullity and finite deficiency) for all values of f3 and 7 except for the values 
— 2 + ~ + k and — 2 — fc where k is any nonnegative integer. In fact, H66|) is an isomorphism for 
— n/q < f3, 7 < — 2 + n/p (recall that we are assuming n > 3, so such /3, 7 exist). Since we are 
principally interested in the behavior of functions at the origin, we will fix 70 e {—n/q, —2 + n/p). 
Then, for (3 e (—2 + n/p, — 1 + n/p), we find that I|tj6|) is surjective with a one-dimensional 
nullspace spanned by |x| 2_ ". 

Now let us consider the formal adjoint of A, which also defines a bounded operator on these 
spaces 

(67) L = n,>.r)i),ih : ^(R"\{0}) i /3+2 , 7+2 (R"), 

where a%j, of course, denotes the complex conjugate of ay. Because a,ij{x) = Sij for \x\ > 1 and 
aij(x) — 5ij vanishes as x — > 0, the analysis in the above references shows that the operator (|67|l 
is Fredholm for exactly the same values of j3 and 7 as for l|tj(j|) . In fact, for fixed nonexceptional 
values of (3 and 7, we may take 8 sufficiently small and use perturbation theory (cf. (T3|, Ch.IV, 
Sec. 5) to conclude that the nullity and deficiency of l|(jtj|) and I|(j7fl agree. 
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So, in addition to the fixed 70 G (—n/q, —2 + n/p), let us now fix (3\ G (—n/q, —2 + n/p) and 
/?2 G (—2 + n/p, —1 + n/p), and denote the adjoints of the corresponding operators H67(l by Ai 
and A2'. 

(68) A 1 : i P _ ft _ 2 ,_ 70 _ 2 (K") - W: 2 p f^ o (R n \{0}) 
is an isomorphism, and 

(69) A 2 : ^ 2 _ 2 ,_ 70 _ 2 (R") -> W:^_ lo (R n \{0}) 

is injective with a one-dimensional cokernel. An arbitrary non-zero functional in Coker^4 2 will 
be denoted by (. 

We introduce a cut-off function 77 G C^°(Bi) equal to 1 on i?i/ 2 . ft follows from l|63f) that 
r]Z G L p _ 0i _ 2 (Bi) but it is not in L p _^_ 2 {Bi). Let F = A{(1 - 7])Z) = -A(r]Z), Since F = 
on B 1/2 and on B\, it follows that F G VFr^ p _ 7o (K™\{0}) C WC^ P _ 7o (»™\{0}). But ® is 
an isomorphism, so nZ is the only solution of A\v = —F. Since L p _p 2 _ 2 (B\) C L V ^_ 2 (B{), 
A2V = —F has no solutions, which means that C(-F) 7^ 0. 

Now let u G Lf oc (Bi) be a weak solution of Au = satisfying M p (u,r) < C r 2 ~~ n+e °. This 
estimate on u implies that u G L p _ l3 _ 2 {Bi) for /3 < —n/q + Eo. So let us restrict our choice of (3\ 
to the interval (—n/q, —n/q + e )- Denote by C an arbitrary constant. Since A(n(u — CZ)) = 
on By 2 and on BJ, we have A(r}(u - CZ)) G Wr|^ 2 _ 7o _ 2 (M n \{0}). Choosing C to satisfy 

can = c(A(nu)), 

we obtain ^(^(77(7/ — CZ))) — 0, which implies 77(7/ — CZ) G Dom(^4 2 ), and in particular 

But this implies that w = u — CZ satisfies M p (w,r) < C rP 2+2 ~~ n / p . Assuming that we had fixed 
/3 2 = j } — 1 — £1 where e± G (0, 1), we conclude that 

(70) M p {w,r) < cr x ~ ei . 
This completes the proof. 
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